Precise Estimation of Cosmological Parameters Using a More Accurate Likelihood 
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The estimation of cosmological parameters from a given data set requires a construction of a 
likelihood function which, in general, has a complicated functional form. We adopt a Gaussian 
copula and constructed a copula likelihood function for the convergence power spectrum from a 
weak lensing survey. We show that the parameter estimation based on the Gaussian likelihood 
erroneously introduces a systematic shift in the confidence region, in particular for a parameter 
of the dark energy equation of state w. Thus, the copula likelihood should be used in future 
cosmological observations. 
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Introduction. Currently, a number of wide-field weak 
lensing (WL) surveys are planned, such as Subaru Hy- 
per Suprime Cam Survey jl] • the Panoramic Survey 
Telescope & Rapid Response System (Pan-STARRS 1 ), 
the Dark Energy Survey (DES 2 ), the Large Synoptic 
Survey Telescope (LSST 3 ) jthe Joint Dark Energy Mis- 
sion (JDEM 4 ), and Euclid [2|, to address questions about 
the nature of dark energy and / or the properties of grav- 
ity on cosmological scales. To achieve the full potential 
of these planned future surveys, it is of great importance 
to employ adequate statistical measures and methods of 
weak lensing for estimating cosmological parameters. 

For the cosmological parameter estimation, almost all 
previous authors have used the \ 2 method in weak lens- 
ing analysis [e.g. [343] ■ However, the probability distribu- 
tion function (PDF) of the weak lensing power spectrum 
was found to be well approximated by a x 2 distribution 
with a heavier positive tail than expected from the nor- 
mal x 2 distribution Q . The % 2 distribution deviates from 
Gaussian distribution on large scales because the num- 
ber of modes corresponding to the degrees of freedom 
are very small. Meanwhile, the % 2 distribution converges 
Gaussian distribution at high I because of the central 
limit theorem. We have to make maximal use of these 
pieces of information accurately to constrain the cosmo- 
logical parameters. If such information is not taken into 
account the likelihood function, the derived cosmological 
parameters can be systematically biased 1(1 U | ■ 

In a companion paper [l2j ]. we constructed a more 



properly accurate likelihood function using the Gaus- 
sian copula (hereafter "copula likelihood") for the cos- 
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mic shear power spectrum, rather than the multivariate 
Gaussian distribution. We show that the copula likeli- 
hood well reproduces the n-dimensional probability dis- 
tribution of the cosmic shear power spectrum estimated 
from 1000 realizations obtained from ray-tracing simula- 
tions performed by 

In this Letter, we estimate the cosmological parame- 
ters using both the copula and Gaussian likelihoods in 
order to evaluate how the difference between the two 
likelihoods affects the parameter estimation. Using the 
Markov Chain Monte Carlo method, we also examine the 
impact on both ongoing and future surveys. 

The cosmological parameters employed for our ray- 
tracing simulations are consistent with the WMAP 3-yr 
results (WMAP3) [HI]. Detailed descriptions of our ray- 
tracing simulations are summarized in Q (see also, (l3|). 

Copulas. The likelihood function plays a central role 
in various statistical analyses. In the companion paper 
[l2| , we derive the copula likelihood for the cosmic shear 
power spectrum. In this section, we briefly summarize 
the method. 

From the Sklar's theorem (l5j . one can relate any n- 
point cumulative probability distribution (CDF) to one- 
point CDFs as 

Prob(xi < xi,x 2 < x 2 , ■■ ■ ,x n < x n ) = F(xi,x 2 , ■■■,x n ) 
= C(F 1 (x 1 ),F 2 {x 2 ),...,F n (x n )). (1) 

Here Xi (i = 1, 2, . . . , n) are independent and identically- 
distributed observed variables (in this work, the shear 
power spectrum, P K (£)), C denotes the function called 
a copula, F denotes the n-point CDF, and Fi denotes 
the one-point CDFs (in this work, the cumulative % 2 dis- 
tributions). Thus, the copula describes how one-point 
CDFs are joined together to give an n-point CDF. A 
comprehensive proof of Sklar's theorem and rigorous def- 
inition of a copula are found in 3- lij ■ 
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The multivariate Gaussian copula is a copula of n- 
dimensional random vector that is multivariate normally 
distributed. This copula is expressed as 

C(u 1 ,u 2 ,...,u n ) ^(^(ifi),^ 1 ^),...,^ 1 ^)) , 

(2) 

where $ x 1 is the inverse function of a one-point Gaussian 
CDF with mean /i^ and standard deviation <ji, and m = 
Fi(xi). Here $ is an n-point Gaussian CDF defined by 



<&(xi,x 2 , ...,x n ) 



— CO J — OO 
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— OO ^(27r)"det(Cov) 



x exp ( —-Ax — /x) T Cov 1 {x — fi) ) dxidx 2 . . .dx, n 



(3) 



with mean fi and n x n covariance matrix. Cov - 
shows the inverse covariance matrix. Hereafter, \x = 
(fii,fi2, ■ ■ ■ ,fin), x = (xi,x 2 , ■ ■ ■ ,x n ), and superscript T' 
stands for the transpose of vector. 

Using this copula, a log- likelihood ln£ is derived as 

m 

n n 

-2 In £ c (xi,X2, ■ • -,x„) = - /U^Cov -1 ^ - fij) 

«=i i=i 

{q% - Hi) 2 



i=l 



2^1n/ l (x l ), 



for a general probability distribution, while for Gaussian 
probability distribution it reduces 

n n 

-2hi£ g (xi,X2, ...,£n) = y^y^ (^i-/ii)Cov~ 1 (f J — /ij). 

(5) 

Here, functions /j represent the derivatives of the CDFs 
Fi, and we have omitted an irrelevant constant term in 
the above two equations. The relation between Ui and qi 
is 



qi = (T l ^ 1 1 (u i ) + fi u 



(6) 



where f i is a cumulative standard normal distribution. 
In what follows, we use Eqs. (@} and Eq. ([5]) to constrain 
the cosmological parameters and investigate how the dif- 
ferent likelihood functions cause a difference in the pos- 
terior distribution function of cosmological parameters. 

Methodology and parameter choices. In the case 
considered in this Letter, the observed variables x i are 
binned values of the nonlinear convergence power spec- 
trum P K (£i), which are computed using the prescription 
of [n| for the nonlinear effects. Note that we assume 
a bin width Aln£ = 0.3, single source redshift distribu- 
tion, i.e. all lensed galaxies lie at z s = 1.0, and we do 



not consider intrinsic ellipticity dispersion er c through- 
out this Letter. Instead, we use the information up to 
multipole £ max = 1071, because the weak lensing power 
spectrum estimated from a realistic survey on smaller 
scales than £ ~ 1000 is expected to be contaminated 
by the intrinsic ellipticity noise [20j . When constructing 
the copula likelihood function, /ij and Cov -1 are esti- 
mated from ray-tracing simulations performed by [9] . It 
is known to be appropriate to adopt a \ 2 distribution for 
the general probability distribution fi, because the one- 
point PDF of the convergence power spectrum is fairly 
well described by % 2 distribution with mean and variance 
of P K (£i) = (P K (eA) and a 2 UA = (P^) 2 ) - P«(4) 2 , 
respectively [see, 9|. From [2l|, this % 2 distribution is 
written as 



p K (i, 



\T-1 



r(T) 



,-p K (ii)/p K {ti) 
PSi) 



(7) 



for P R (li) > and fa = for P K {li) < 0. Here T(x) 
is the gamma function and we define T = P K {ti) 2 j <J 2 (li) 
which corresponds to the number of independent modes. 

For simplicity, we work with two cosmological param- 
eters in weak lensing given as 



P = (^dr 



,lnl0 10 A2 



(8) 



where ttdm is the dark matter density today and A R is the 
amplitude of primordial curvature perturbations defined 
(4) at k = 0.002 Mpc 1 whose fiducial values in our 



analysis are 0.196 and 20.45 x 10~ 10 , respectively. The 
other cosmological parameters are fixed at the fiducial 
values of the WMAP 3-yr cosmology. When we estimate 
the equation of state parameter w in the dark energy 
equation, combined with another probe of WMAP3 data, 
we work with three cosmological parameters given as 



p= (r» dm , W ,lnl0 10 A2). 



(9) 



The parameter space we explore is as follows: fidm = 
[0.002,1.0], w = [-6.0,0], lnl0 10 A| = [-2.5,8.0] and 
the fiducial values of these parameters are taken as 0.196, 
— 1 and 3.02, respectively. Assuming flat priors for the 
cosmological parameters, we employ the Markov Chain 
Monte Carlo method [22| to constrain cosmological pa- 
rameters given the cosmological observables. Sixteen 
parallel chains were computed and the convergence test 
is made based on the Gclman and Rubin statistics called 
tl R — 1" statistics 23j. In this work, each chain typically 
has 300,000 points and R - 1 < 0.02 for both of the two 
models. 

Results. In this section, we constrain the cosmologi- 
cal parameters using both the copula and Gaussian like- 
lihoods and evaluate how the difference between the two 
likelihoods affects the estimation. We consider two cases: 
ongoing and future surveys. 
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Impact on ongoing weak lensing surveys. In our ray- 
tracing simulations, the survey area is set as fi s = 25 
deg 2 . Therefore the fundamental mode of ray-tracing 
simulations is if = 72. 

Figure Q] shows la and 2a confidence level constraints 
on the fl m -as plane. Note that fl m denotes the present- 
day total matter density, i.e., the sum of the current dark 
matter and baryon densities. The red and blue contours 
show the marginalized constraints obtained by the cop- 
ula likelihood (Eq. 2]) and Gaussian likelihood (Eq. [5]), 
respectively. Clearly the results from the copula likeli- 
hood (red contours) and Gaussian likelihood (blue con- 
tours) are significantly different. The contours obtained 
with the copula likelihood are shifted toward a parame- 
ter region that gives a lower convergence power compared 
to that from the Gaussian likelihood. This result is at- 
tributed to the fact that the median of the x 2 distribution 
for the convergence power spectrum is smaller than that 
of the Gaussian distribution. 

Figure [2] shows the constraints on the fl m -w plane ob- 
tained when the weak lensing information is combined 
with data from WMAP3. Also in this case we see the 
difference between the results from the two likelihoods. 
Similarly, the contours obtained by the copula likelihood 
are shifted toward the parameter region that gives the 
lower power, i.e. smaller Q m and larger w, than the 
Gaussian likelihood, as expected. Consequently, if one 
uses the (approximate) multivariate Gaussian likelihood, 
one will overestimate the lower bound and underesti- 
mate the upper one for the constraint on w. Specif- 
ically, we found an allowed range of the equation of 
state parameter w for the Gaussian likelihood model as 
— 1.055 < w < —0.956 , while for the copula likelihood 
model as —1.038 < w < —0.941 at 95% confidence level. 
Therefore, a few percent of systematic error potentially 
exists in the parameter estimation if one uses approxi- 
mate Gaussian likelihood. 

Impact on future weak lensing surveys. Next let us 
assume a survey area of 2000 deg 2 , which roughly cor- 
responds to the planned Subaru Weak Lensing Survey 
[lj]. Because our ray-tracing simulation is limited to 25 
deg 2 , we have to extend it to estimate the mean power 
spectrum /ij and their covariance with help from ana- 
lytic treatments. When estimating the mean power spec- 
trum fa, we use the linear perturbation theory based 
on WMAP 3-yr fiducial cosmology for multipoles in the 
range 8 < I < if — 72, where if corresponds to the 
largest angular mode of our ray-tracing simulation (5x 
5 deg 2 for an area). For the covariance matrix Cov, 
we assume that it has only diagonal elements for mul- 
tipoles smaller than if — 72 because the linear approxi- 
mation is valid for that multipole range. In the absence of 
shape noise, the diagonal elements would arise from sam- 
ple variance and should be equal to the power spectrum 
squared divided by the number of independent modes 
in the bin. Also, we assume that the covariances be- 
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FIG. 1: Two-dimensional constraints (la and 2a C.L.) on 
Cl m and as for survey area S7 S = 25deg 2 . The red and blue 
contours show the constraints obtained by likelihood based on 
the copula likelihood (Eq.[4]) and Gaussian likelihood (Eq. [5]), 
respectively. The cross symbol shows the fiducial cosmological 
parameters. 
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FIG. 2: Same as Fig. [T] but two-dimensional marginalized 
constraints on f2 m and w for weak lensing plus WMAP 3-yr 
results. 



tween the scales larger than if = 72 and scales smaller 
than if = 72 are zero, which means that the powers at 
the largest scales are independent from those at smaller 
scales. 

Figure |3] shows the constraints on the (ft m , w) plane 
obtained when weak lensing information with the sur- 
vey area of 2000 deg 2 is combined with WMAP3. In 
this case, we do not see any significant difference be- 
tween the results from the two likelihoods. This result 
is understood as follows. If one considers the larger sur- 
vey area, T becomes larger at a fixed multipole because 
a 2 becomes smaller (see, Eq. [7]). Hence, \ 2 distribu- 
tion for the convergence power spectrum becomes well 
approximated by the Gaussian distribution at the fixed 
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FIG. 3: Same as Fig. [2] but assumed survey area is 2000 deg 2 
for weak lensing. This good agreement may not be caused by 
a physical reason, but rather attributed to the coarse binning. 

multipole. Even though at the lowest multipole bin the 
distribution may be deviated from Gaussian, the bulk of 
information comes from the large multipole bins, where 
the distribution is almost Gaussian. Hence the contours 
become indistinguishable. 

It should be noted, however, that the Gaussianity in 
the distribution of the power spectrum should be due 
to our sparse binning of the multipole space. To utilize 
the information as much as possible, we have to divide 
the multipole space into a larger number of bins than we 
have done in this study. In such a case, the distribution 
of the power spectrum at each multipole bin will deviate 
from the Gaussian distribution and the copula likelihood 
will still be useful to model the probability of the power 
spectrum distribution. 

Summary. In this Letter we constructed a likelihood 
function for the convergence power spectrum from weak 
lensing survey using copula. The copula likelihood func- 
tion describes more accurately the distribution of the 
power spectrum derived from ray-tracing simulations and 
thus utilizes more information than the previous meth- 
ods in the literature. We then used the copula likelihood 
model to calculate the allowed range of the cosmological 
parameters paying particular attention to the dark en- 
ergy equation of state parameter w, and compare it with 
that derived from the conventional multivariate Gaussian 
likelihood model. We found that, for the 25 deg 2 weak 
lensing survey, the results can be different even combined 
with the CMB data, depending on which likelihood func- 
tion is used. The difference is as large as a few percent. 

For the 2000 deg 2 weak lensing survey, we found the 
results coincide with each other. This result allows us 
to use a multivariate Gaussian likelihood model for the 
future weak lensing survey, which will greatly simplify 
the parameter estimation analysis. We note, however, 



that this coincidence might be an artifact from our sparse 
binning of the multipole space. A full analysis with finer 
binning is difficult at present because computational cost 
is high, and we leave it for a future study. 
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